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ABSTRACT 

We show that recently published assertions that advection dominated accretion 
flows (ADAFs) require the presence of strong winds are unfounded because they as- 
sume that low radiative efficiency in flows accreting at low rates onto black holes im- 
plies vanishing radial energy and angular momentum fluxes through the flow (which 
is also formulated in terms of the 'Bernoulli function' being positive). This, however, 
is a property only of self-similar solutions which are an inadequate representation of 
global accretion flows. We recall general properties of accretion flows onto black holes 
and show that such, necessarily transonic, flows may have either positive or negative 
Bernoulli function depending on the flow viscosity. Flows with low viscosities (a 0.1 
in the a-viscosity model) have a negative Bernoulli function. Without exception, all 
2-D and 1-D numerical models of low viscosity flows constructed to date experience 
no significant outflows. At high viscosities the presence of outflows depends on the 
assumed viscosity, the equation of state and on the outer boundary condition. The 
positive sign of the Bernoulli function invoked in this context is irrelevant to the pres- 
ence of outflows. As an illustration, we recall 2-D numerical models with moderate 
viscosity that have positive values of the Bernoulli function and experience no out- 
flows. ADAFs, therefore, do not differ from this point of view from thin Keplerian 
discs: they may have, but they do not have to have strong winds. 
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1 INTRODUCTION 

In many systems containing accreting galactic and extra- 
galactic black holes, the luminosity deduced from observa- 
tions is much lower than the one obtained by assuming a 
'standard' radiation efficiency of ~ 0.1. It has been pro- 
posed that accretion in such systems can be modeled by 
advectively dominated accretion flows - ADAFs (for recent 
reviews see Kato, Fukue & Mineshige 1998; Abramowicz, 
Bjornsson & Pringle 1998; Lasota 1999). In ADAFs, the 
main cooling process is advection of heat — radiative cool- 
ing is only a small perturbation in the energy balance and 
has no dynamical importance. ADAFs are quite successful 
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in explaining observed spectral properties of accretion onto 
black hole in low mass X-ray binaries, Galactic center, and 
some active galactic nuclei. 

Recently, Blandford & Begelman (1999, hereafter 
BB99) put in question the physical self-consistency of ADAF 
models by arguing that flows with small radiative efficiency 
should experience outflows, which in some cases could be 
so strong as to prevent accretion of matter onto the black 
hole. Low luminosities, therefore, would not be due to low 
radiative efficiency, but simply to absence of accreting mat- 
ter near the central object. Hence, ADAFs should be re- 
placed by what BB99 call ADIOs, i.e. 'advection-dominated 
inflows-outflows' (later 'advective' has been changed to 'adi- 
abatic'). The BB99 argument about the necessity of outflows 
from ADAFs is based on two statements about what they 
call the 'Bernoulli constant' (the significance of this 'con- 
stant' - in reality a function - is discussed in the next Sec- 
tion): [1] The Bernoulli constant in flows with low radiative 
efficiency must always be positive, a point which was first 
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made by Narayan & Yi (1994), [2] A positive Bernoulli con- 
stant implies outflows. 

In the present paper we show that [1] the vertically 
integrated Bernoulli function is everywhere negative in small 
viscosity ADAFs which have a vanishing viscous torque at 
the flow inner 'edge', and are matched to the standard thin 
Shakura-Sunyaev disc (SSD) at the outer edge, [2] A positive 
Bernoulli function in an ADAF does not imply outflows — 
in a representative class of 2-D numerical models of ADAFs 
with moderate viscosity, Bernoulli function may be positive 
but despite of that outflows are always absent. 

The main conclusion of our paper is that while it is not 
yet clear whether some ADAFs with high viscosity could in- 
deed have significant outflows, both general theoretical ar- 
guments and numerical simulations point out that it is un- 
likely that ADAFs with low viscosity could experience even 
moderate outflows of purely hydrodynamical origin. Similar 
conclusion has been recently obtained by Nakamura (1998), 
who used a different approach. 

We use here cylindrical coordinates (r, z, (f>) and denote 
gravitational radius of the accreting black hole (with the 
mass M) by rG = 2GM/c 2 . We model the gravitational field 
of the black hole by the Paczyhski & Wiita (1980) potential. 
We assume in this paper that viscosity is small. This means 
that a 0.1 for the kinematic viscosity coefficient given by 
a phenomenological formula 



Cti pC 'sound, 



(1.1) 



where lp — P/|VP| is the pressure length scale and C SO und 
is the sound speed. Our theoretical arguments do not de- 
pend on a particular functional form (prescription) of a. 
In numerical models we assume the 'standard' prescription, 
a = const. Recent theoretical estimates based on numerical 
simulations of turbulent viscosity (see e.g. Balbus, Hawley 
& Stone 1996) show that most likely ct 0.1. In this pa- 
per we shall consider mainly this low viscosity range of a. 
Arguments based on fitting predicted spectral properties of 
ADAF models to observations (see Narayan 1999 for review) 
seem to require the moderate viscosity range, 0.1^q$ 0.3, 
that partially overlaps with the low viscosity range consid- 
ered in this paper. As pointed out by Rees (private commu- 
nication) it is not clear if the high viscosity range, a <; 0.3, 
is physically realistic. 



2 THE BERNOULLI CONSTANT, FUNCTION 
AND PARAMETER 

In stationary, inviscid flows with no energy sources or losses, 
the quantity, 



1 9 

B = W + -V + $ = const, 



(2.1) 



is constant along each individual stream line, but, in general, 
is different for different stream lines. This quantity is called 
the Bernoulli constant. Here W is the specific enthalpy, V 
is the velocity (all three components included), and <& is the 
gravitational potential. Obviously, a particular streamline 
may end up at infinity only if Bo > along it. The exis- 
tence of stream lines with Bo > is therefore a necessary 
condition for outflows in stationary inviscid flows with no 
energy sources or losses, and Bo < for all streamlines is 
a sufficient condition for the absence of outflows. However, 



Bo > is not a sufficient condition for outflows. For ex- 
ample, in the case of the classical Bondi's accretion (pure 
inflow) Bo is a universal positive constant. 

In all viscous flows, Bo defined by (2.1) is, of course, 
not constant along individual stream lines. However, the so- 
called 'Bernoulli parameter', 



1 K 



w + iv 2 + $ 



(2.2) 



is a universal constant in 1-D, vertically integrated, self- 
similar Newtonian models of ADAFs introduced by Narayan 
& Yi (1994). Here Vk is the Keplerian velocity. The reason 
for that is simple: in self-similar models all quantities scale 
as some power of the cylindrical coordinate r. In particular, 
because both 



B = W + iy 2 + $ / const, 



(2.3) 



and Vj| scale as r _1 , their ratio -Bo must be a constant. 
Narayan & Yi (1994) and BB99 have argued that because 
the vertically integrated models have Bo > 0, in generic 2-D 
flows there should be streamlines escaping to infinity. 

This argument is not correct for the simple reason 
that neither Bo > (as already mentioned), nor Bo > 0, 
are sufficient conditions for outflow existence. Note that 
Bo = constant is not a real physical conservation law, but 
rather an artifact induced by a purely mathematical assump- 
tion that the flow is self-similar. This very strong assump- 
tion is motivated only by practical convenience, not by the 
physical properties of a flow. Indeed, in all numerical, global 
models of ADAFs constructed so far, self-similarity does not 
represent well the global flow properties and Bo changes 
its value and even its sign from place to place. Large scale 
changes are due to the global balance between viscous heat- 
ing, advective cooling and PdV work. Small scale changes 
are typical for ADAFs in which a strong convection (circula- 
tion) is present (these ADAFs have a moderate or small vis- 
cosity) . In this case, the sign of Bo changes rather abruptly, 
tracing convective bubbles. In addition, at a given place, the 
value and the sign of Bo strongly fluctuates in time. 

The Bernoulli constant Bo is not a useful quantity for 
ADAF study, because they are viscous flows. The Bernoulli 
parameter Bo has no physical significance and therefore it 
is not a convenient quantity for theoretical arguments. We 
shall use here the 'Bernoulli function', which is formally de- 
fined by (2.3), but — of course — in ADAFs (or any other 
viscous flows) it is not constant along streamlines. As in the 
case of steady dissipation-free flows, B < everywhere is a 
sufficient condition for the absence of outflows. 

We were forced to introduce such new object 217 years 
after the death of Daniel Bernoulli because of the above- 
mentioned articles on outflows from ADAFs, which use the 
related concepts of the Bernoulli constant or 'parameter'. It 
is easier, as shown below, to refute their arguments by using 
the same 'paradigm'. 



3 THE INNER BOUNDARY CONDITION 

In this Section we would like to stress the importance of the 
fact that the accreting body is a black hole, in particular, the 
implications of the transonic nature of the accretion flow. 
We will assume that the flow viscosity is low. Paczynski 
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(1978, unpublished) noticed that in this case the mass loss 
from the inner part of an accretion disc around a black hole 
is fully governed by the general relativistic effect of 'rela- 
tivistic Roche lobe overflow' which operates close to a 'cusp' 
radius r cusp , where the angular momentum of the disc takes 
the Keplerian value, j(r cuap ) = jic(rcusp) (see Abramowicz 
1981, 1985 for details). Here jk = Vkt is the Keplerian an- 
gular momentum. The cusp is formed by a critical equipo- 
tential surface, similar (but topologically different) to the 
critical Roche surface in the binary stellar-system problem. 

Jaroszyhski, Abramowicz & Paczyhski (1980) proved 
that for small viscosity accretion discs one must have r m b < 
Tcuap < r ms , where r m t, and r ma are respectively the radii 
of the marginally bound, and the marginally stable Kep- 
lerian circular orbits around a black hole. In the case of 
a non-rotating (Schwarzschild) black hole, as well as in 
the Paczyhski & Wiita (1980) potential, r m (, = 2rc and 
r m , = 3r<3. At r ms the Keplerian angular momentum has 



and 



a minimum value j„ 



(3/2) 



3/2 



r G c. 



Numerical simula- 



tions show that low viscosity accretion discs with reason- 
able outer boundary conditions (discussed in the next Sec- 
tion) have always a high specific angular momentum, that 



is for r c 



< r < ro, they are always super-Keplerian: 



j(r) > ji<(r), where r ~ 5rc (see Figure 1). 

The mass loss through the cusp, at the rate M cusp , 
induces an advective cooling Qc Uap which is a very sensi- 
tive function of the vertical thickness h of the disc at the 
cusp: independent of the equation of state (i.e. indepen- 
dent of the adiabatic index of the accreted matter 7) one 
has Qcus P ~ £/i 3 , where E is the disc surface density at 
the cusp. This implies that the relativistic Roche lobe over- 
flow stabilizes possible unstable thermal modes in the region 
close to r CU s P because an overheating would cause vertical 
expansion, which then would induce strong advective cool- 
ing. More precisely, in terms of Pringle's thermal stability 
criterion (Pringle 1976; Piran 1978), 



dlnh 



In Q cusp 
dlnh 



(3.1) 



where Q^ is ~ T,h 2 is the rate of viscous heating, one has 
2 = (l.h.s) < (r.h.s) = 3, which proves thermal stability 
(Abramowicz, 1981). Thus, an increase in the Bernoulli func- 
tion caused by overheating does not produce outflows, but 
enhances the inflow into the black hole. This analytic pre- 
diction has been confirmed in all details by 2-D time depen- 
dent numerical simulations. In particular, Igumenshchev, 
Chen & Abramowicz (1996) showed that the analytic for- 
mula (Abramowicz, 1985) for the rate of the mass inflow 
into the black hole induced by the relativistic Roche lobe 
overflow reproduces, in a wide range of parameters and with 
impressive quantitative accuracy, the behavior of M CUBp cal- 
culated in all their numerical simulations. 

For low viscosity (high angular momentum) black hole 
accretion flows, Abramowicz & Zurek (1981) found that the 
regularity condition at the sonic point V 2 = C 2 ound requires 



C sou 



£ 2 3k(T sonic) j (r sonic) ^ ^ 



r 2 ■ C 2 

some 



Because e < 1 and j K (r 'sonic 

j {t sonic) 



c c, one has 



JK (V sonic} 



= 1 



e 2 + 0(e 4 ), (3.3) 



= 1-0(0, 



(3.4) 



i.e. the sonic point is very close to the cusp. These properties 
have been confirmed by 1-D and 2-D numerical simulations 
performed independently by numerous authors. 

The supersonic flow at r < r SO nic does not hydrody- 
namically influence the subsonic flow at r > r sonic- For this 
reason, somewhere at an 'inner edge' r in ~ r son4c ~ r cusp , 
the viscous torque vanishes, g(ri n ) = 0. From (2.3) and 
(3.2) - (3.4) one derives, for a polytropic fluid, i.e. with 
W = C 2 ound /{ 1 -l), 



4 ( % S O Ti 



+ 



1)2 2(7-1) 



e 2 + 0(e 4 ). (3.5) 



Here x = r/ra- In the physically relevant region, 2 < 
Xsonic < 3, the leading term (zeroth order in e 2 ) of this 
function equals, independently of 7, to negative binding en- 
ergy at circular Keplerian orbit and varies between and 
-1/16. 

Thus, for ADAFs with low viscosity (high angular mo- 
mentum) one should adopt the following inner boundary 
condition, 



g{r in ) = 0, j(r in ) = j K {r ln ). 



(3.6) 



because this is imposed by fundamental properties of the 
black hole gravitational field and thus must be always 
obeyed. In addition, if at the sonic point C 'sound /c 1, the 
sonic point regularity condition demands that the Bernoulli 
function at the inner edge should be negative, 

B{r ln ) w B(r S omc) < 0, (3.7) 

independent of the equation of state (independent of 7). 

4 THE OUTER BOUNDARY CONDITION 

ADAFs cannot exist for arbitrary large radii. For example, 
according to Abramowicz et al. (1995) ADAFs cannot extent 
beyond the radius 



C^-r G > Tin, 



(4.1) 



3k i T sonic) 



where m = M/MEdd and C ~ 10 2 . In general, C and the 
power of a depend on the cooling mechanisms included into 
the model (see Menou, Narayan & Lasota 1999). 

Observations suggest that, in several systems the in- 
ner ADAF is surrounded by a geometrically thin, standard 
Shakura-Sunyaev disc so that there must exist a transition 
region where, for r ~ r ou t the ADAF is matched to a Keple- 
rian disc. The physical mechanism which triggers the transi- 
tion is not known (see, however, Meyer & Meyer-Hofmeister 
1994; Kato & Nakamura 1998) and it is not clear what is 
the relation between r max and r out . 

These uncertainties make the conditions at the outer 
edge r ou t of low viscosity accretion discs less precisely de- 
termined than the conditions at the inner edge. Despite of 
this, Abramowicz, Igumenshchev & Lasota (1998) found a 
simple analytic argument that proves that independent of 
the physical reason for the transition, the angular momen- 
tum close at r ou t should have exactly the Keplerian value, 
j{r ut) ~ 3K{r ou t). This was confirmed by numerical models 
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Figure 1. Angular momentum distribution for low viscosity 
ADAF (solid line) that fulfills inner (3.6) and outer (4.2) bound- 
ary conditions, given by Paczynski's fitting formula (5.2). The 
dashed line corresponds to the Keplerian distribution j}((r). 



A-ou 



f(r)-f K (r) 



dr 



Ji j%(r)-f(r) 



dr, 



and thus the 1/r weighted area A ollt between these curves 
in the region [ro,r ut] is also roughly determined. In addi- 
tion, one should have dj/dr > and d(j /r 2 )/dr < 0. 

The function j(r) can be approximated by the analytic 
fitting formula used in the Paczynski's (1998) toy model for 
ADAFs, 



j'M = jri. 



1+6 



\ Tin / \ Tin / 



1.5/q 



(5.2) 

Note that this formula is different from the one in the pub- 
lished version of Paczynski's article. The error it contained 
has been corrected in the astro-ph version. Here a and 6 are 
constants that depend on r; n and r out , and can be deter- 
mined from (5.1). Although in actual models the range of 
radii for which the angular momentum is approximately con- 
stant is much reduced, one should stress, that (5.2) is a good 
qualitative representation of the generic angular momentum 
distribution in ADAFs when the magnitude of viscosity is 
low independent of the functional form of viscosity. 



of the only specific ADAF-SSD transition model worked out 
to date, in which the transition occurs due to the presence 
of a turbulent flux (Honma 1996, Kato & Nakamura 1998). 
With an ADAF joining the SSD at the outer edge, one 

has, 

j(r out ) = jK(r ut), B(r out ) ~ e K {r out ) < 0, (4.2) 

Here eK{r out ) is the negative Keplerian orbital binding en- 
ergy. The second equation in (4.2) follows from 

B(r out ) = e K (r out ) [l + O ((h out /r out ) 2 )] < 0, (4.3) 

where h ou t <C r out is the thickness of the outer SSD. In all 
numerical simulations of low viscosity, high angular momen- 
tum accretion flows, the outer boundary condition (4.2) is 
approximately fulfilled at all radii r 3> r in independent of 
whether the ADAF-SSD transition occurs. 



5 THE ANGULAR MOMENTUM 

The shape of angular momentum distribution between r in 
and r ut depends mainly on viscosity but, as numerous 1-D 
and 2-D models demonstrate, in the low viscosity case it is 
always similar to that shown in Figure 1. The physical rea- 
son for such a shape is clear. Close to the inner edge, the 
viscous torque is ineffective (note that g(ri n ) ~ 0) and thus 
the specific angular momentum has a very small gradient. 
Thus, the location of r^ n roughly determines the location 
of the first crossing point ro between j(r) and jk (r) curves, 
and therefore also the 1 /r 3 weighted area Ai n between these 
curves in the region [n n , ro] . Mechanical equilibrium condi- 
tion demands that (Abramowicz, Calvani & Nobili 1980), 



6 THE BERNOULLI FUNCTION 

Imagine a cylinder r = const crossing an ADAF from its 
upper, z = h(r), to lower z = —h(r), surface. Let Mo, Jo and 
Eo denote, respectively, the total amount of mass (Mo < 0), 
angular momentum and energy that cross the surface of the 
cylinder per unit time. In a steady state, from the Navier- 
Stockes equations of mass, angular momentum and energy 
conservations integrated along the cylinder it follows that, 

r+h(r) 

Mo = / 2-Krp(r, z)V r (r, z)dz = const, (6.1a) 

J -h(r) 

Jo = Moj(r) +g{r) = const, (6.16) 

E = M B(r) + (l(r)g(r) = const. (6.1c) 

Here fi(r) = j(r)/r 2 is the angular velocity. The radiative 
energy flux from ADAFs is very small and it was ignored 
in (6.1c). Derivation of (6.1) assumes that the flow has an 
azimuthal symmetry (no dependence on <f)) , and that the or- 
bital velocity Vj, is much greater than the accretion velocity 
V r (and 'vertical' velocities) which is true, except very close 
to the inner edge, for flows with a small viscosity considered 
in this Section. Except h(r), each radial function that appear 
in (6.1) represents the averaged value of the corresponding 
quantity. In particular, 

B(r) = — ■ / 27vrp(r,z)V r (r,z)B(r,z)dz. (6.2) 

2h(r)M J_ h(r) 

The same equations (6.1), with the same assumptions and 
with the same averaging procedure (6.2), have been used by 
BB99, and by Paczyhski (1998) in his recent illuminating 
paper on a toy model of ADAFs. 
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According to Narayan and Yi (1994) in self-similar mod- 
els of ADAFs, both j(r) and g(r) scale as r 1 ^ 2 , and both 
B(r) and Q,(r)g(r) scale as r~ . These scaling properties 
imply that j — Cjr 1 ^ 2 and E — Cet -1 , with Cj = const, 
Ce = const. Thus, the fluxes Jo and Eo can be constant if 
and only if C.j = Ce = 0. This implies that both the angular 
momentum flux, and the energy flux are exactly zero in self- 
similar models: Jo = 0, and Eo — (BB99). From the first 
of these equations it follows that g(r) = — j(r)/Mo- This, 
together with the second equation show that the Bernoulli 
function must be positive, 

B(r) = n(r)j(r) = [rQ(r)] 2 > 0. (6.3) 

The conclusion B(r) > is based on self-similar solu- 
tions which are obviously inconsistent with boundary con- 
ditions. We shall show that by taking properly into account 
the boundary conditions (3.6) and (4.2), one arrives at the 
opposite conclusion: B(r) < 0. 

The inner boundary condition g(ri n ) — implies that 

Moj(r in ) = M j(r) + g(r), (6.4a) 

MoB(r in ) = M B(r) + fi(r)<?(r). (6.46) 

From the last two equations one derives 

B{r) = B{r m ) + V{r)j{r) - Q(r)j(r in ). (6.5) 

Equation (6.5) yields, at the outer edge of the disc r out 

B(r in ) = B(r out ) - Sl{r ou t)\j{r ut) - j(r i7l )]. (6.6) 

Because in stable discs j(r out ) > j(ri n ), the last term on the 
right hand side of this equation is always negative. From this 
fact and from equation (4.3) one concludes that for standard 
ADAFs, i.e. for those that have a vanishing torque at the 
inner edge ti n and match the standard Shakura-Sunyaev 
disc at the outer edge r out one must have 

B(r ln ) < B(r out ) < 0. (6.7) 

Thus, the Bernoulli function in standard ADAFs must be 
negative both at the inner and outer edges. Identical con- 
clusions have been reached, but not explicitly stated, by 
Paczyhski (1998): see his equation (20) which is equivalent 
to our equation (6.6). 

From equation (5.2) describing a typical angular mo- 
mentum distribution in the disc, and equations (6.5), (6.6) 
one may calculate B(r) in the whole disc. Figure 2 shows the 
Bernoulli function by the thick solid line. The thin broken 
line shows the prediction of the self-similar model. In the 
same Figure 2 we present for comparison by the thin solid 
line the time-averaged Bernoulli function for 2-D numerical 
model of ADAF with a = 0.01. Details of numerical tech- 
nics described by Igumenshchev & Abramowicz (1999). The 
model has r in — 3ro, r out = 8000rc, and 7 = 5/3. 

One concludes that in ADAFs with small viscosity, 
which fulfill standard boundary conditions, the Bernoulli 
function must be everywhere negative. This conclusion is 
independent of the functional form of viscosity. Obviously, 
a flow which has the Bernoulli function that is everywhere 
negative does not experience outflows. Note, however that 
B(r) calculated here is averaged with respect to the flux of 
mass. Thus, it may happen that B > in some polar di- 
rections. This is indeed the case close to the disc surface of 
some low viscosity numerical models calculated by ICA96. 
These models show weak outflows, with M out <C Mq. 
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Figure 2. The Bernoulli function in ADAFs with small viscosity 
is everywhere negative. The solid thick line corresponds to the 
angular momentum distribution given in Figure 1, with r;„ = 
2.003rG, r ou t = lO^rg. The thin dashed line corresponds to the 
Narayan & Yi (1994) self-similar solution (with 7 = 3/2), not 
compatible with the boundary conditions, used by BB99. The 
thin solid line repersents the time-averaged Bernoulli function for 
2-D numerical model of ADAF with a = 0.01. 

7 NUMERICAL SIMULATIONS 

Analytic arguments presented in the previous Section and 
pointing out that no significant outflows of hydrodynamical 
origin should be present in low viscosity ADAFs have been 
fully confirmed by all numerical simulations performed to 
date. Below we give a list of some representative works. 
1-D global simulations of transonic flows in optically 
thick case (slim discs): Abramowicz, Czerny, Lasota & 
Szuszkiewicz (1988); Kato, Honma & Matsumoto (1988ab); 
Chen & Taam (1993); Szuszkiewicz & Miller (1997). 

1- D global simulations of transonic flows in the optically 
thin case (ADAFs): Honma (1996); Chen, Abramowicz & 
Lasota (1997); Narayan, Kato & Honma (1997); Gammie 
& Popham (1998); Nakamura, Kusunose, Matsumoto & 
Kato (1998); Igumenshchev, Abramowicz & Novikov (1998); 
Ogilvie (1999). 

2- D time dependent simulations of ADAFs: Igumenshchev, 
Chen & Abramowicz (1996); Igumenshchev & Abramowicz 
(1999); Stone, Pringle & Begelman (1999). 



8 ADAFS WITH LARGE VISCOSITY 

We have seen that there is a significant theoretical and nu- 
merical evidence which shows that purely hydrodynamical 
outflows in ADAFs with low viscosity are unlikely to occur. 
Physical reasons for the absence of outflows in this case seem 
to be well understood. They depend on some fundamental 
properties of the black hole gravity. 

The situation in the case of moderate and large viscosity 
is less clear. The analytic calculation of B(r) presented in 
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the Section 5 cannot be repeated in the case of large viscosity 
because in this calculation one assumes that V r <C V^>, while 
in flows with large viscosity, all velocity components are of 
the same order. This brings an additional unknown term to 
the energy equation (6.1c) so that the system has too many 
unknowns to be solved. Also, one can not further assume 
zero viscous torque acting at Vi„ , if non-circular motions are 
significant. Thus, equations (6.1) are insufficient to calculate 
B(r) even in the self-similar case. In the large viscosity case 
one can argue neither that self-similarity implies B(r) > 0, 
nor that the boundary conditions imply B(r) < 0. 

Purely hydrodynamic outflows have been seen in nu- 
merical 2-D simulations of high viscous accretion flows (Igu- 
menshchev & Abramowicz 1999), and in 2-D self-similar 
models (Xu & Chen 1997, and references therein), but they 
are not a universal property of viscous flows. From the exist- 
ing results it is obvious that the presence of outflows depends 
on the magnitude of viscosity parameter a, adiabatic index 
7, and, probably, on the outer boundary conditions at r ou t. 
However, the boundary J-{a, 7, ...) that divide the parame- 
ter space into the flows with and without outflows is yet to 
be found. 

Certainly, one cannot argue that B > implies out- 
flows. We illustrate this point by showing in Figure 3 a model 
of ADAF (a = 0.3, 7 = 3/2, r in = 3r G , r out = 8000r G ) in 
which B > and no outflows. The model has been calcu- 
lated using numerical technics described by Igumenshchev 
& Abramowicz (1999). In Figure 4, for the same model, we 
show calculated distributions of B(r) (solid line), 



B(r) = 



27T7^ 

M c 2 



GM \ 



normalized viscous energy flux (dashed line), 



G(r) = 



"i^ 2 Jo 



(VrUrr + V 9 U r g + V^Rr*) COS 8d6, 

(8.2) 



and normalized total energy flux (solid thick line), 
E (r) 



M c 2 



= B(r) + G(r). 



(8.3) 



In (8.1)-(8.3) r, 9 and <j> are spherical coordinates and II is 
the shear stress tensor. In the model the inward energy ad- 
vection flux [term B in (8.3) and short-dashed line in Fig. 4] 
almost equals to the outward viscous energy flux [term G 
in (8.3) and long-dashed line in Fig. 4] at each radius. This 
behaviour is similar to that in the self-similar Narayan & 
Yi (1994) solution, where the oppositely directed fluxes ex- 
actly compensate. For comparison, we show B(r) for the 
self-similar solution in Figure 4 by dotted line. When cor- 
rect boundary conditions are taken into account there is no 
exact compensation and the total energy flux Eq must be a 
small (in absolute value) non-positive constant in the sta- 
tionary dissipative accretion flow of the type discussed here. 
Note that in the actual model Eq (solid line in Fig. 4) is not 
constant and oscillates with a small amplitude due to an 
inaccuracy of our numerical scheme, which does not exactly 
conserve the total energy balance. This inaccuracy is inside 
« 5% of relative error, which is acceptable for our purposes. 




Figure 3. 2-D numerical stationary ADAF model with a = 0.3, 
7 = 3/2, ri n = 3tq and r ou t = 8000rG- Only inner part of the 
model is shown. Arrows indicate the magnitude of the vector r 2 pv, 
where r is the spherical radius. Regions with positive Bernoulli 
function are in white, while regions with negative values are in 
grey. Clearly, B > is not sufficient for outflows. 



\ 

\ 

\ 

\ 

\ 

. \ 






\ \ 

\ \ 






- f ~~ — 






/ 

/ 

- / 
/ 






/ 

- / 

- / 






10 


100 


1000 




r/r r 





Figure 4. Radial profiles of 0-averaged quantities, given by eqs 
(8.1)-(8.3), for the model presented in Figure 3. Profiles of B(r) 
(short-dashed line), G(r) (long-dashed line) and Eo(r)/ Mqc 2 
(solid line) are shown. The dotted line corresponds to B(r) for 
the self-similar Narayan & Yi (1994) solution. See text for expla- 
nations. 
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Absence of winds in ADAFs 7 



9 CONCLUSIONS 

1. Significant outflows of purely hydrodynamical origin are 
not present in ADAFs with low (a ^ 0.1) viscosity. This 
conclusion follow from general theoretical arguments that 
involve fundamental properties of black hole gravity, and 
are well understood. All numerical simulations confirm this 
theoretical arguments and show no outflows. 

2. It would be very interesting to perform a systematic in- 
vestigation in the parameter space {a, 7, r out } by filling it 
with models of ADAFs and find the regions (necessarily with 
a > 0.1) corresponding to outflows. 

3. Outflows from ADAFs might occur due to non- 
hydrodynamical factors such as magnetic fields, radiation, 
etc. These processes cannot be modeled in purely hydrody- 
namical terms by adopting special value of 7 or form of the 
Bernoulli function: the extra physics should enter through 
solutions of the relevant equations (see e.g. King & Begel- 
man 1999). 
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